Using Absorption Imaging to Study Ion Dynamics in an Ultracold Neutral Plasma 
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We report optical absorption imaging of ultracold neutral plasmas. Images are used to measure the 
ion absorption spectrum, which is Doppler-broadened. Through the spectral width, we monitor ion 
equilibration in the first 250 ns after plasma formation. The equilibration leaves ions on the border 
between the weakly coupled gaseous and strongly coupled liquid states. On a longer timescale of 
microseconds, we observe radial acceleration of ions resulting from pressure exerted by the trapped 
electron gas. 



Plasma physics traditionally studies systems with tem- 
peratures of thousands of kelvin or greater because col- 
lisional ionization of atoms requires kinetic energies on 
this scale. Ultracold neutral plasmas, created by pho- 
toionizing laser-cooled and trapped atoms, access an ex- 
otic regime in which particle energies can be on the order 
of 1 K. 

Fundamental interest in ultracold neutral plasmas 
stems from a range of phenomena in the ultracold regime. 
Recent experiments have studied plasma creation [J, col- 
lective modes |2|, and recombination into highly excited 
Rydberg atomic states @- The spontaneous evolution 
of a dense, ultracold gas of Rydberg atoms into an ul- 
tracold plasma has also been investigated p|. This is 
part of a search in atomic clouds for an analog of the 
Mott insulator-conductor transition in condensed mat- 
ter [jj. Recombination in these systems resembles the 
methods used to produce cold antihydrogen with trapped 
positrons and antiprotons 

A series of theory papers 0,0, E3, El [13 explored is- 
sues surrounding thermalization and recombination in ul- 
tracold neutral plasmas when both electrons and ions are 
near or in the strongly coupled regime In strongly 
coupled plasmas the electrical interaction energy between 
the charged particles exceeds the average kinetic energy. 
This reverses the traditional energy hierarchy that un- 
derlies our normal understanding of plasmas based on 
concepts such as Debye screening and hydrodynamics. 
Strongly coupled plasmas exist in dense astrophysical 
systems |l4j . matter irradiated with intense laser fields 
\wi Il6l| . dusty plasmas of highly charged macroscopic 
particles 0], an d non-neutral trapped ion plasmas [3 
that are laser cooled until they freeze into Wigner crys- 
tals. 

We now report the first results with a new probe of 
ultracold plasmas: absorption imaging. This technique 
provides in situ, non destructive measurements, and of- 
fers excellent spatial, temporal, and spectral resolution. 
We describe the use of this probe to study ion-ion equili- 
bration and expansion of the plasma during the first few 
microseconds after photoionization, but we emphasize its 
great potential to study a host of phenomena such as ion 
collective modes ^^|, shock waves 0], recombination, 
and particle-particle spatial correlations |18j . 

The production of an ultracold neutral plasma starts 
with atoms that are cooled and confined in a magneto- 
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FIG. 1: Atomic and ionic energy levels involved in the exper- 
iment, with decay rates. (A) Neutral atoms are laser cooled 
and trapped in a magneto-optical trap (MOT) operating on 
the 1 So — 1 Pi transition at 461 nm. Atoms excited to the *Pi 
level by the MOT lasers are ionized by photons from a laser 
at ~ 412 nm. (B) Ions are imaged using the 2 Si/2 — 2 P\ 
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transition at 422 nm. P 1 / 2 ions decay to the -D3/2 state 7% 
of the time, after which they cease to interact with the imag- 
ing beam. This does not complicate the experiment because 
ions typically scatter fewer than one photon during the time 
the imaging beam is on. 



optical trap (MOT) (Fi g. 1 A). This aspect of the exper- 
iment was described in j2f)j | . The neutral atom cloud is 
characterized by a temperature of about 20 mK and a 
density distribution given by n(r) — n exp(— r 2 /2a 2 ), 
with cr=l mm and n — (4 ± 2) x 10 10 cm" 3 . The 
number of trapped atoms is (6 ± 1) x 10 8 . 

To form the plasma, the MOT magnets are turned off 
and atoms are ionized with photons from a 10 ns pulsed 
dye laser whose wavelength is tuned just above the ioniza- 
tion continuum (Fig. 1 A) . Because of the small electron- 
to-ion mass ratio, the initial electron kinetic energy (E e ) 
approximately equals the difference between the photon 
energy and the ionization potential. E e /kB can be as 
low as the bandwidth of the ionizing laser, which is ^100 
mK. The initial kinetic energy for the resulting singly- 
charged, electronic ground state ions is close to that of 
the original neutral atoms. As we will discuss below, the 
resulting non-equilibrium plasma evolves rapidly. Up to 
12 ± 1% of the neutral atoms are ionized producing plas- 
mas with a peak density of (5 ± 3) x 10 9 cm" 3 . 

Immediately after photoionization, the charge distri- 
bution is neutral everywhere. Due to the kinetic en- 
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FIG. 2: Optical depth of an ultracold neutral plasma. (A) 
The delay between the formation of the plasma and image 
exposure is 100 ns, and the initial peak density is not = 
5 x 10 9 cm -3 . (B) Residuals of a fit to a two-dimensional 
Gaussian profile are close to the photon shot-noise. 

ergy of the electrons, the electron cloud expands on 
the timescale of the inverse electron plasma frequency 
T e = cjpg 1 = y/ m e eo/n e e 2 < 1 ns, where m e , n e , and 
e are the electron mass, density and charge. On this 
timescale the ions are essentially immobile. The result- 
ing local charge imbalance creates a Coulomb potential 
energy well that traps all but a small fraction (< 5%) 
of the electrons. Simulations 0] show that electrons es- 
cape mostly from the edges of the spatial distribution, 
and the center of the cloud is well described as a neutral 
plasma plj . The diagnostic used in previous experiments 
was charged particle detection of electrons and ions after 
they had left the plasma. 

Spectroscopic diagnostics are ubiquitous in plasma 
experiments, and some even provide spatial informa- 
tion jBuch as spatially resolved laser induced fluorescence 
[lSl |22| ■ The absorption imaging reported here is par- 
ticularly well-adapted for small, cold, relatively dilute 
plasmas that evolve very quickly. It is also a powerful 
techni que for studying laser cooled and trapped neutral 
atoms j23j. A collimated laser beam, tuned near reso- 
nance with the principle transition in the ions (Fig. IB), 
illuminates the plasma and falls on an image intensified 
CCD camera. Following Beer's law, the optical depth 
{OD) is defined in terms of the image intensity without 
({background) and with (Ipiasma) the plasma present, 

OD(x,y) = \n(I background (x,y)/I p i 

/CO 
-00 

where nni is the peak ion density, and a(v) is the absorp- 
tion cross section at the image beam frequency, v. In the 
last line we have inserted a Gaussian density distribution 
for the ions, which leads to the function used to fit the 
data. 

Figure [21 shows a typical absorption image. The inten- 
sity of the probe beam is approximately 200 /iW/cm 2 , 
which is much less than the saturation intensity of the 
transition (38 mW/cm 2 ). The spatial resolution of typ- 
ical images is about 100 /im, limited by pixel averaging 



performed to improve the signal-to-noise ratio of the im- 
ages. 

To study the time-evolution of the plasma, we vary 
the delay between the formation of the plasma and im- 
age exposure (t delay) with 10 ns accuracy. The minimum 
camera exposure gate width is 50 ns. For the shortest 
exposure times, which we use for the best time resolu- 
tion at very short tdeiay, we typically average about 60 
ionizing laser shots. For longer delay times we use longer 
image exposure times of up to 800 ns, and decrease the 
number of accumulations to keep the total number of 
photons detected approximately constant. The repeti- 
tion rate for ionization and image recording is about 5 
times per second. 

The imaging beam itself can also be gated with a min- 
imum width of about 200 ns. We use this capability for 
tdeiay > 1 /xs to turn off the image beam until the camera 
exposure begins. This prevents optical pumping of the 
ions to the 2 D 3 / 2 state (Fig. IB). 

Plotting the peak optical depth as a function of im- 
age laser frequency provides the absorption spectrum of 
the ions (Fig. |3J) . The imaging laser linewidth of about 
5 MHz is negligible on the scale of the 21.5 MHz natural 
linewidth of the transition. As described below, addi- 
tional broadening of the absorption spectrum provides a 
wealth of information on the plasma dynamics. 

Here we describe experiments using the time evolution 
of the absorption spectrum to study ion dynamics for a 
plasma with Ni = 7 X 10 7 ions, initial peak density for 
ions and electrons of n oi as n 0e — (5 ± 3) x 10 9 cm~ 3 , 
and E e /kB — (68 ± 5) K. We chose a relatively large E e 
in order to avoid complications that arise when the elec- 
tron Coulomb coupling parameter (r e = e 2 /47T£n aksTe) 
approaches or initially exceeds unity, such as screening 
of the ion interaction |24|, and rapid collisional recom- 
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FIG. 3: Absorption spectra of ultracold neutral plasmas. We 
plot the peak optical depth derived from fits to data such 
as Fig. 2. The frequency is with respect to a Doppler-free 
absorption feature in a strontium discharge cell. Both spectra 
correspond to E^/ks = 68 K and the same initial peak plasma 
density of not = 5 x 10 9 cm -3 . Data are fit with Voigt profiles, 
and the increase in linewidth for longer tdeiay is clear. 
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FIG. 4: Ion dynamics determined from the spectral width. 
Error bars are two sigma uncertainties from Voigt profile fits 
to absorption spectra. The time after photoionization corre- 
sponds to the timing of the center of the camera exposure 
gate. Data show disorder-induced heating of the ions as the 
ions equilibrate during tdeiay < 250 ns. (Inset) For longer 
times, the Doppler broadening of the spectrum continues to 
increase as ions accelerate radially due to electron pressure. 
This process is described by a hydrodynamic model with no 
free parameters. 



bination and heating of the electrons 0, 0, 0, Il2| . 
Here, a = (47rn 0e /3) _1 / 3 is the Wigner-Seitz radius, 
and T e — ^E e /k B is the electron temperature as set by 
the wavelength of the ionizing laser. For this sample 

r e = o.i. 

The observed spectral linewidths are all significantly 
broader than the natural linewidth of the transition. Col- 
lisional broadening is not significant because of the rela- 
tively low particle density and collision frequency in the 
sample. The dominant contribution to the linewidth is 
Doppler broadening, which makes the spectrum a very 
accurate probe of the ion velocity distribution. From the 
fit of each spectrum to a Voigt profile, we extract the 

where to; is 



rms Doppler broadening, <td = B ^ 
the ion mass, and A is the wavelength of the transition. 
This allows us to determine an ion temperature (Ti) as a 
function of time, as shown in Fig. 

The temperature increases rapidly for tdeiay < 250 ns. 
Two pieces of information imply that this is thermal- 
ization of ions with themselves after creation in a spa- 
tially disordered state. The timescale is on the order 
of the inverse plasma frequency of ions, Tj = cj" 1 = 

\J miEQ / riQie 2 = 100 ns, which is the timescale on which 
ions respond to perturbations from their equilibrium spa- 
tial distribution. The temperature of 1.4 ± 0.1 K derived 
from the Doppler width at 250 ns is also on the order 
of the amount of potential energy inherent in the initial 
disorder of the ions. A simple estimate implies there is 
e 2 /AnsoakB ~ 5 K of potential energy that will be redis- 
tributed during thermalization. 

This thermalization has been modelled with molecu- 



lar dynamics simulations 0, 0, ^| , and the experimental 
data is in general agreement with the theory. The accu- 
racy of the imaging probe will enable detailed comparison 
of measurement and theory for the time dependence of 
the ion temperature. 

The final temperature reached, however, can be com- 
pared with an expression derived in Assuming com- 
plete initial disorder and incorporating the screening ef- 
fects of the electrons, 
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Here, k = 
(e a k B T e /n 0c 
quantity U = 



a/X D 

,2)1/2 
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= V3r e = 0.55 where Ad 
7 jjxa is the Debye length. The 
is the potential energy per parti- 



cle in units of e 2 / Altera. It has been studied with molec- 
ular dynamics simulations |24| for a homogeneous sys- 
tem of particles interacting through a Yukawa potential, 
(j){r) = 4 ^ exp(— r/Aj), which describes ions in the 
background of weakly coupled electrons j2^| . The exper- 
imental values for rii, Ti, and k, imply U—-0.73. Equa- 
tion |2 then gives 7j = 1.4 K, in excellent agreement with 
the measured value. 

We now address the level of Coulomb coupling for the 
ions. From the temperature and peak density we derive 
r.; = 3±1 for the thermalized ion cloud at tdeiay — 250 ns. 
For a system of charges embedded in a uniform neutraliz- 
ing background, formally called a one-component plasma 
[26j , local spatial correlations characteristic of a strongly 
coupled fluid appear for T > 2. For Debye shielded parti- 
cles, such as the ions studied here, screening will reduce 
the correlations. This effect is approximately incorpo- 
rated by using the effective coupling constant T* = Te~ K 
d Q. At 250 ns, r* = 2 ± 1, and the ions are just 
on the edge of the strongly coupled fluid phase. Perhaps 
other initial exp erimental parameters, or laser cooling of 
the ions @, l27j . will lead to more strongly coupled sys- 
tems, although r* = 2 already puts the experiment in an 
interesting regime. 

For tdeiay > 250 ns, the spectral width continues to 
increase, but at a slower rate, as shown in the inset of 
Fig. 0] This slow increase is due to outward radial accel- 
eration of the ions caused by pressure exerted by the gas 
of trapped electrons. This was studied experimentally in 
and theoretically by a variety of means in |n.J. The 
experiments measured the final velocity that the ions ac- 
quired, which was approximately yj E e /mi. Here we ob- 
serve the ion dynamics at much earlier times during the 
acceleration phase. 

A hydrodynamic treatment predicts that the force 
per ion is 
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V(n e (r)fc B T e ) _ „rk B T e 



rii(r) 



(3) 



where the ion and electron density distributions are 
n e (r) w tii(r) — noiexp(— r 2 /2cr 2 ). We approximate the 
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rms size (<7j) as the mean of the observed sizes a x and 
Oy , and we assume thermal equilibrium for the electrons 
throughout the cloud 0]. This force leads to a radial 
expansion velocity for the ions, v r (r) , which is correlated 
with position and increases linearly with time. This does 
not represent an increase in the random thermal veloc- 
ity spread or temperature of the ions. Due to the large 
mass difference, thermalization of ions and electrons |9j 
is slow and occurs on a millisecond timescale. The in- 
crease in Doppler broadening due to thermalization is 
approximately one order of magnitude smaller than what 
is observed. 

Using the dynamics implied by Eq. [31 the evolution 
of the Doppler broadening can be calculated. The mean 
squared velocity component along the imaging laser is 

(vl) = J d 3 rdv T ^-P{v T )(vT+v r (r)co S 6) 2 ,(4) 

where P(vt) is the thermal distribution of v z for Ti = 
1. 4 K . The resulting theoretical Doppler broadening 
VY^D/A, plotted in Fig. 0] reproduces the data accu- 
rately. 

The ion acceleration is thus an excellent diagnostic of 
the electron temperature. This will be of great value 
in future studies because the temperature is predicted 
to evolve in a complicated fashion for lower initial T e 
due to recombination and disorder- induced heating of the 
electrons 0, El 0, E3- One also expects that there 
will be a cooling effect at longer tdeiay due to expansion 
of the plasma and evaporative cooling The small 



discrepancy between theory and data for long tdeiay in 
Fig. 01 may indicate the onset of this cooling, although 
expansion of the plasma is small on the timescale of these 
observations. For the maximum tdeiay and typical v r rs 
15 m/s, v r tdeiay ~ 50 /im. This small increase in size is 
observed in the images (Fig. EJ . 

The initial study using absorption imaging of an ul- 
tracold neutral plasma has probed ion dynamics in the 
first few microseconds after photoionization. It revealed 
disorder- induced heating that was predicted in and 
showed that the ions equilibrate on the boundary of the 
gas-liquid transition. Acceleration of ions due to electron 
pressure was also evident, and can be used to monitor the 
electron temperature. 

Many future experiments suggest themselves. Some 
of the most interesting are investigating dynamics when 
the initial electron Coulomb coupling parameter is large 
and recombination and disorder-induced electron heating 
are expected to dominate the plasma evolution. Detailed 
study of ion and electron thermalization at the border 
of the strongly coupled regime is also possible. Improve- 
ments in the imaging optics will significantly increase the 
image signal-to-noise ratio and allow study of features on 
the ion density distribution with ~ 10 /jm experimental 
resolution. 
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